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Most general relativity textbooks devote considerable space to the simplest example of a black 
hole containing a singularity, the Schwarzschild geometry. However only a few discuss the dynamical 
process of gravitational collapse, by which black holes and singularities form. We present here two 
types of analytic models for this process, which we believe are the simplest available; the first involves 
collapsing spherical shells of light, analyzed mainly in Eddington-Finkelstein coordinates; the second 
involves collapsing spheres filled with a perfect fluid, analyzed mainly in Painleve-Gullstrand coordi- 
nates. Our main goal is pedagogical simplicity and algebraic completeness, but we also present some 
results that we believe are new, such as the collapse of a light shell in Kruskal-Szekeres coordinates. 



PACS numbers: ll.10.Kk, 11.25.Uv, 11.30.Pb, 98.80.Es 

I. INTRODUCTION 

Black holes and singularities are certainly some of 
the most peculiar and interesting features of general 
relativity[l|. From the perspective of fundamental the- 
ory their properties are far from completely understood, 
especially their relation to quantum theory 2]; from the 
perspective of astrophysics and observational astronomy 
they seem more and more to play a central role in the 
physics of objects ranging from stellar to cosmological 
size 3]. 

Most general relativity textbooks devote considerable 
space to discussing the simplest example of a black hole 
containinga singulari ty, the Schwarzschild geometry^ 
H H H H H H [lfl [H O , and some also discuss Kerr's 
generalization to a rotating system, although only a few 
derive it 0,0]. Most also explain, at least qualitatively, 
why black holes should exist in the real world, based 
on the instability of neutron stars of more than about 
4 solar masses; even neutron degeneracy pressure is not 
sufficient to provide stability. A few texts also discuss 
Schwarzschild's interior solution for a constant density 
starUQ, which clearly illustrates some of the same in- 
stability features as more realistic neutron star models. 
But the actual dynamical process of collapse, whereby 
a massive body becomes a black hole, is a more com- 
plex dynamical problem, and is often either neglected or 
treated heuristically and qualitatively. 

For a brief history of gravitational collapse see 
reference fli). The study of gravitational collapse dynam- 
ics began with the classic seminal work of Oppenheimer 
and Snyder [l^ , who set up the equations and gave a 
semi-quantitative discussion of the collapse of a spherical 
non-radiating star. They also found a complete analytic 
solution for the case of a uniform perfect fluid of zero 
pressure - now widely referred to as "dust." This exam- 
ple is less artificial than it might appear since the effect 
of pressure in a more realistic treatment turns out to be 
not very critical; but the treatment by Oppenheimer and 
Snyder involves some awkward algebra in its treatment 



of the boundary between the interior and exterior of the 
dust star. 

Some textbooks discuss collapse much like Oppen- 
heimer and Snyder, notably Landau and Lifshitz[f| and 
Misner Thorne and Wheeler 0. In particular these texts 
use a Novikov type coordinate system, in which the 
dust is co-moving in a bound orbit |16||. This coordinate 
system is synchronous and co- moving, and is wonder- 
fully simple conceptually, but almost equally awful alge- 
braically. Due to the algebraic difficulty, we believe, this 
is not the simplest way to treat the problem. Most other 
texts give a heuristic but convincing analysis by compar- 
ing the motion of a dust particle just outside a collapsing 
dust sphere with the motion of the surface; matching this 
motion is equivalent to matching the exterior and inte- 
rior geometries, which is a rather obvious but important 
theorem [l7j. 

What we have tried to do in the present paper is to 
analyze dynamical collapse in as simple a way as pos- 
sible by our choice of coordinates; we use Eddington- 
Finkelstein [H [3 

type coordinates for infalling incoher- 
ent light or "null dust," and Painleve-Gullstrand^ 
type coordinates for perfect fluids, including dust. These 
coordinates seem almost magically well suited to this 
purpose, as we hope will become apparent in sections 
2 and 5. We use zero energy fluid systems in contrast to 
negative energy bound systems as used in references [3] 
and [5], largely because the algebra is much simpler 
This choice seems to us well justified because the zero 
and nonzero energy systems have the same qualitative 
behavior as they approach the black hole state as seen 
from outside in terms of standard Schwarzschild time, 
and have exactly the same behavior as they approach 
the singularity as seen in terms of the proper time of a 
falling observer. Moreover the zero energy case does not 
have the awkward conceptual problem of the motion be- 
fore the maximum size is reached, since the maximum 
size occurs at an infinite time in the past. We there- 
fore believe the zero energy case is more illustrative and 
economical of effort. 
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Almost needless to say we have treated only spheri- 
cally symmetric collapsing systems. The Birkhoff theo- 
rem guarantees that a spherically symmetric system does 
not emit gravitational radiation; non-sphcrically sym- 
metric collapse generally entails gravitational radiation, 
making the problem vastly more complex, though also 
more interesting. 

To achieve our goal of simple textbook type examples 
of collapse, we must pay a price, which is some degree of 
artificiality. The main artificiality is the spherical sym- 
metry; also the in-falling shells of light or other massless 
material of sections 2 and 3 are not likely to be found 
in nature; similarly we do not expect to find the perfect 
fluid sphere surrounded by an elastic shell of section 5, 
or the zero pressure perfect fluid of sections 6 and 7. 

Since our purpose in the present paper is overtly ped- 
agogical, much of what we obtain here is known and 
available somewhere in the vast research literature, al- 
though obtained with different techniques and different 
coordinates. However we believe our treatment is the 
simplest, both conceptually and algebraically, mainly due 
to our coordinate choice. We do not know of any use of 
Painleve-Gullstrand coordinates in the present context, 
although they are becoming more widely used[22L |23| : 
as will be seen in section 6 these coordinates combine a 
Schwarzschild type radial coordinate and a Friedmann- 
Robertson- Walker (FRW) type cosmological time coordi- 
nate to give a natural interpolating system for describing 
both the interior and exterior of a fluid sphere. While 
the Eddington-Finkelstein and Painleve-Gullstrand coor- 
dinates serve their mathematical purpose quite well, they 
have a number of drawbacks such as a "time" marker 
that is not well-behaved everywhere in spacetime, and 
they are not synchronous; we have therefore also treated 
the collapse of a thin shell of light by transforming to 
Kruskal-Szekeres coordinates, which are conformally flat 
and thus display the causal structure of the collapse pro- 
cess quite clearly |3 E3- In Kruskal-Szekeres coordi- 
nates one has a nice picture for distinguishing real world 
black holes from "eternal" black holes, with their associ- 
ated wormholes and white hole segments, none of which 
may be expected to actually exist in our universe. 

We use throughout this paper a rather novel technique 
for handling the boundary conditions between the fluid 
and the exterior, which we believe is physically clear. 

Since the pioneering work of Oppcnhcimer and Snyder 
there has been an enormous amount of work on gravi- 
tational collapse and related topics, most directed to re- 
search rather than pedagogy. The present pedagogical 
paper originated from several research problems, one in- 
volving a formation mechanism for singularity-free black 
holes filled with heavy vacuum, and another involving a 
finite version of the standard or "concordance" cosmolog- 
ical model. For these we needed analytic collapse models 
that were algebraically simpler and more explicit than 
we found in the literature. Our hope is that our methods 
and models might also be of use to students and teachers 
in illustrating more standard gravitational collapse. 



This paper is organized as follows: the first part (sec- 
tions 2 to 4) deals with collapsing shells of incoherent 
light or null dust, which is entirely characterized by 
its null 4-velocity and energy density; the second part 
(sections 6 and 7) deals with spheres of perfect fluid, 
including some with nonzero pressure and some with 
non-uniform density. In section 2 we use Eddington- 
Finkelstein coordinates to obtain the metric for the col- 
lapse of a thin light shell onto a black hole to form a larger 
black hole, and the collapse of a thin light shell to form 
a black hole ab initio. The ab initio collapse is amusing 
in that it is almost certainly the simplest complete sce- 
nario for the formation of a black hole, albeit it a rather 
artificial one. In section 3 we use the results of section 
2 to obtain the metric for the collapse of a thick shell of 
light, layer by layer. The density profile of the shell is 
largely arbitrary, and we give one specific example. In 
section 4 we transform the metric for the thin light shell 
collapse to conformably flat Kruskal-Szekeres coordinates 
in order to display the causal structure of the geometry 
in a clear way - nearly equivalent to a Penrose diagram; 
as already noted we use these coordinates to emphasize 
the difference between black holes formed by collapse and 
eternal Schwarzschild black holes. In section 5 we present 
our way of doing the classic problem of collapse of a uni- 
form fluid sphere, using Painleve-Gullstrand coordinates 
for maximum simplicity. We also obtain an explicit solu- 
tion for a nonzero pressure fluid with a linear equation of 
state: p = ap\ to balance the pressure gradient force at 
the surface we use the artifice of a thin elastic shell with 
tangential pressure instead of the slowly decreasing radial 
pressure of a real stellar system; our model is essentially a 
balloon with uniform internal density and pressure. The 
stress energy tensor of the surface is described by the ra- 
dius of the sphere, its energy density, and the parameter 
a; for a = we recover the standard results for the dust 
ball. In section 6 we deal with zero pressure non-uniform 
dust spheres, with the density characterized by a largely 
arbitrary function; again using Painleve-Gullstrand co- 
ordinates we build up the system layer by layer, to ob- 
tain the standard results of Oppenheimer and Snyder and 
Landau and Lifshitz etc. We do not allow the dust layers 
to cross, which could giv e rise to infinite densities and su- 
perficial singularities 26] . We give one specific example 
of the metric of a non-uniform sphere; the special case 
of uniform density agrees with the analysis of section 5. 
Since this paper is mainly pedagogical we have through- 
out sacrificed brevity and included considerable algebraic 
detail; only elementary mathematical methods are used. 

It is worth mentioning that some simple and fairly ob- 
vious extensions of the present work can be made. First, 
by a reversal of time the collapse of light shells can be 
viewed as outgoing radiation, so we may easily obtain 
results like those of Vaidya for the geometry of a radiat- 
ing body [22j]. Similarly the collapse of the perfect fluid 
can be time reversed to yield the metric for a black hole 
emitting matter, that is a type of white hole. Also it is 
straightforward to include a cosmological constant term 
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in the equations to obtain the metric for a collapsing sys- 
tem in an exterior Schwarzschild de Sitter geometry; due 
to its somewhat more cumbersome algebra we have not 
included this in the present paper, and leave it as an exer- 
cise. (See reference 28].) We will say more about further 
extensions and applications in section 8. 



II. THIN LIGHT SHELLS 

The motion of light in Schwarzschild geometry is alge- 
braically awkward at and inside the Schwarzschild radius 
when Schwarzschild coordinates are used, but it becomes 
quite simple in Eddington-Finkelstein (EF) coordinates. 
In this section we use EF coordinates to study the mo- 
tion of light, the growth of a black hole by absorption of 
light, and the ab initio formation of a black hole by a thin 
shell of light. The last process is most likely the simplest 
model for the formation of a black hole by gravitational 
collapse. 

The Schwarzschild metric in Schwarzschild coordinates 
is, with c = 1, 



ds 2 



(1 + u)dt 2 s 



dr 2 
1 + 



l dtt 2 



(1) 



where u = -R/r and d£l 2 = d9 2 + sin 2 9d(j) 2 . The black 
hole surface is a sphere at the Schwarzschild radius (twice 
the geometric mass) R = 2m = 2GM . It is both an infi- 
nite redshift surface where 500 = and a null surface or 
one-way-membrane. In terms of the Schwarzschild time 
coordinate both light and particles take an infinite time 
to reach the black hole surface from the exterior. Specif- 
ically, for light falling radially from 



At s = n 



Rln 



( Tj - R \ 
\r-RJ 



(2) 



Thus light never reaches the surface but approaches 
asymptotically with a characteristic time R. It is for 
this reason that we introduce the EF time coordinate, in 
terms of which the surface is reached in a finite time. The 
EF time coordinate is obtained from the Schwarzschild 
time by a radially dependent shift 

t s = t + g{r). (3) 

This leads to the EF form of metric, 

ds 2 = dt 2 -dr 2 ~r 2 dfl 2 + u(dt±dr) 2 , (4) 

provided that we choose the transformation function g to 
obey 



dg = u 
dr 1 + u 



The solution to this is 



(•5) 



(6) 



which has an infinite stretch at r — R. Henceforth we 
refer to the coordinates and metric form in Eq.(@J) as EF 
for any function u(r, t). 

The EF form is extremely convenient because radially 
infalling light behaves quite simply if the plus sign in 
Eq.(£l} is chosen. For such light we set the line element 
equal to zero, to obtain 

ds 2 = dt 2 - dr 2 + u(dt + dr) 2 

= (dt + dr) (dt -dr + udt + udr) = 0. (7) 



so that 



dt + dr = 0, t + r = const, 
dr/dt = (1 + u)/(l - u), 



infalling light, (8) 
outgoing light. (9) 



Thus the path of infalling light is independent of the met- 
ric function u, and is the same as in flat space. See Fig.^ 
Conversely, outgoing light "stalls" where u = — 1 and 
.9oo = 0, that is at the infinite redshift surface. 





FIG. 1: In (a) a thin light shell falls into a black hole to produce 
a larger black hole. In (b) a thin light shell with flat Minkowski 
interior collapses to form a black hole. 

We will soon need the Einstein tensor for the EF met- 
ric form Ea.l|Hl. Its nonzero components are easily cal- 
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Here a dot denotes a time derivative and a prime denotes 
a radial derivative. We define the stress energy tensor in 
terms of the Einstein tensor by using the field equations 



a = — 07T(jrl p. 



(ID 



For pure Schwarzschild geometry the Einstein tensor is 
everywhere zero, as is easily verified. 

We now study the formation of a black hole by a sin- 
gle thin shell of light, and in the next section we will 
consider a thick shell of light. For both purposes we be- 
gin with a thin shell of light falling into a Schwarzschild 
black hole with radius R to form a black hole with radius 
R 1 as shown in Fig. the light shell obeys r + t = 
everywhere. We may write the metric in all of spacetime 
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using a step function and its complement 9 = 1 — G , 
as 



R R' 

u = e(r + t) Q(r + t), 



(12) 



which is of course time dependent. Substituting this into 
Ea. (|55|l and Eq.l^TJ we obtain the stress-energy tensor 
for the thin light shell, which is singular due to the step 
function, 



(13) 



where k a = (1,-1,0,0) and hp = (1,1,0,0). The null 
vector k a corresponds to infalling light. From the energy 
density Tq we may calculate the energy or effective mass 
of the light shell by going to a time in the distant past 
when the shell was in asymptotically flat space, 

M s = j ^r 2 T%dr = = M' - M. (14) 

This verifies that the mass of the initial black hole plus 
the energy or effective mass of the light shell equals the 
mass of the final black hole. 

The above results hold for the special case when the 
initial black hole is replaced by flat space, or R = 0. Note 
that everything above is consistent with the interior of a 
spherical thin light shell being flat Minkowski space. This 
represents the formation of a black hole by the gravita- 
tional collapse of a single thin light shell (or other mass- 
less material), as shown in figure lb; this would seem to 
be the simplest complete example of gravitational col- 
lapse. 



III. THICK LIGHT SHELLS 

Our results from the preceding section may be used to 
construct a model for the collapse of a thick light shell, 
layer by layer. The initial state is a sequence of concen- 
tric thin shells with a flat Minkowski interior, as shown 
in Fig.^. Each region of spacetime has a Schwarzschild 
geometry as discussed in the previous section. The inner- 
most shell with energy Ami collapses to form a black hole 
of mass mi = Ami, followed by others to form interme- 
diate black holes of mass m 7 - = Ami + 



+ Arrij , ending 
with a final black hole with mass m/ and Schwarzschild 
radius R = 2m/. Between the shells the geometry is 
given by Eq.Q, with a metric function 



-2mj/r. 



(15) 



The evolving infinite rcdshift surface, defined by u = — 1, 
is the zigzag line in Fig. 

The continuous analog of the sequence is a family of 
concentric shells with a continuous label A, which we 
choose to run from to 1. The parameter A plays the 
same role as the discrete label j: that is the energy or 



effective mass inside the A shell is denoted by m(A). As 
in Ea. l|15[) the metric function in the light shell is 



u = — 2m(A)/r. 



(16) 



The evolving infinite redshift surface is the smooth line 
in Fig. 




k = 1, r+ t = A 
infinite redshift, 
u= -1 



*.=0, r+t = A 



FIG. 2: In (a) a discrete sequence of light shells forms a black hole; 
(b) shows the continuous version of the same process. 

A labeling scheme, which is convenient for both the 
present light shell and the fluid systems to be discussed 
later, is to take the energy inside the A shell to be pro- 
portional to A, 



m(A) = Am/, m/ = final mass. 



(17) 



The time at which the A shell reaches the center may be 
chosen almost arbitrarily as a function of A, subject to 
the constraints that it be for the innermost shell A = 
and increase monotonically to the final time A for the 
outermost shell A = 1, as shown in Fig. [5] We denote 
this function as F _1 (A)A, for reasons that will become 
apparent. The equation of motion for the shell A is thus 



r + t = F-\X)A. 



(18) 



From its definition F 1 must have an inverse, so we may 
invert this to obtain 

A = F (-T-)' m(X) = m f X = m f F, (19) 

with F(0) = and F(l) — 1. Thus F serves as a density 
profile function, and the metric in the interior of the light 
shell may be written as 



2 m/ fr + t 



RF 



(20) 



The infinite redshift surface, defined by u = — 1, is thus 
determined by 



r = RF 



r + t 



(21) 



Inverting this we obtain for the infinite redshift surface 

t = AF-^r/R) - r. (22) 

This implies that if the total duration of collapse A is 
sufficiently large then t must be a positive monotonic 
function of r. 
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Using the metric function Eq. (|20() we may calculate the 
Einstein tensor and the stress-energy tensor from Eq. 1|1(J|) 
and Ea. lfTT|l . This gives for the light shell, 



RF' 



87rGr 2 A 



k a kp, light shell, 



(23) 



where F 1 denotes the derivative of F with respect to its 
argument, and k a is the null vector defined in Ea. (|f 31 . 
From this we obtain the energy density and the total 
energy of the light shell as 



RF' 



8irGr 2 A' 



n-T° - 
P — 1 o — 

o 

M s = / ATtr 2 pdr = 



m f F (!) _ "V f0A s 
G ~ G ' [ ' 



where the integral is again done in the asymptotically flat 
spacetime of the distant past, verifying that the energy 
of the light shell is equal to the final black hole mass. 

Our approach to boundaries and boundary conditions 
is somewhat unorthodox; we do not impose boundary 
conditions per se on the metric function between regions 
of spacetime. Instead we use the above solutions (in vac- 
uum and within the light shell) to calculate the stress- 
energy tensor, which is defined by the field equations and 
Ea. (|10|l . The resulting stress energy tensor must be zero 
in vacuum, correctly represent the stress energy within 
the light shell, and describe the stress on the boundaries. 
If it does, the solution makes physical sense. For the 
region near the inner and outer boundaries of the light 
shell we may write the metric function as 



-(RF/r)G(r + t), 



inner surface, 



-(RF/r)@(r + t— A), outer surface. 



(25) 



Calculating the stress energy tensor from Ea. (|lt)(l and 
Eci. (|llfl we find that there is no singular shell at these 
boundaries since the relevant singular derivatives cancel; 
thus the stress-energy tensor associated with the bound- 
aries is zero and the solution is thus physically reason- 
able. 

As a specific example let us take the density profile 
function to be linear, 



F 



r + t 



m(A) = m f(^-) 



(26) 



This corresponds to a constant rate of energy impact- 
ing the center. From Eq.J3U|) the metric function in the 
various regions is, 

0, Minkowski region, 

// = < —R/r, Schwarzschild region, (27) 

-R(r + t)/(rA), within light shell. 

The infinite redshift surface, from Eq. l|21|l . is the linear 
function 



(28) 



The slope of this is positive for A/R > 1, that is when 
the light shell thickness is greater than its Schwarzschild 
radius. Within the light shell the energy density is, from 
Eq^ and Ea.lfTTJl. 



T° - n - 
1 o — P — 



R 



8nGr 2 A : 



light shell interior, 



(29) 



which is independent of time. 

In Fig. |3 we show a qualitative sketch of some outgoing 
light rays for a fairly general light shell. The sketch is 
made by noting that, from Eq.@: the slope of outgoing 
rays is 1 in the Minkowski region and at large distances 
it approaches 1; the slope is along the infinite redshift 
surface; finally, the slope is -1 at the center of the black 
hole. There is a last-ray-out emitted from the origin to 
infinity, after which all outgoing rays (as well as particles) 
are trapped within the surface at R and eventually fall 
into the singularity. The surface defined by the last-ray- 
out is thus a global horizon. This illustrates that the 
horizon and the infinite redshift surface are quite different 
inside the time dependent light shell, unlike the situation 
for the time independent Schwarzschild geometry. 




FIG. 3: Some "outgoing" light rays in the thick light shell collapse. 
The last ray out (B) hovers at the Schwarzschild radius and defines 
a horizon. 

As a further application of our methods we note that 
our results can be reversed in time to describe radiation 
being emitted by a spherically symmetric system such 
as a star or black hole [H H| S E3 To do this we 
use the minus sign for the metric in Eq.(4) and turn the 
diagram in Fig. 2b upside down. It is then easy to modify 
our algebraic results to describe any reasonable light shell 
density profile. This may be useful in studying the final 
stages of black hole evaporation, when the gravitational 
field of the radiation becomes comparable to that of the 
black hole and cannot be neglected [22^: it might help 
in determining if a black hole radiates entirely away to 
vacuum or leaves behind a remnant 1321 l33l . 



IV. THIN LIGHT SHELLS IN 
KRUSKAL-SZEKERES COORDINATES 

The thin light shell discussed in section 2 is probably 
the simplest model of gravitational collapse. However 
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in EF coordinates goo is negative for the interior of the 
black hole. This means that if one is at coordinate rest 
(dr = dO = d<p = 0) then the square of the proper time 
interval ds 2 has the opposite sign of the square of the 
coordinate time interval dt 2 , so that t may not be inter- 
preted as a good time marker in that region of spacctime. 
The same is true of Schwarzschild coordinates. Kruskal- 
Szekeres (KS) coordinates, which are discussed in many 
texts, were developed to solve this problem^, 0, EU- 
Here we will obtain KS type coordinates for the thin light 
shell collapse of section 2. 

By KS coordinates we mean a system in which the t, r 
part of the metric is conformal to flat Minkowski space, 
with no singularities or zeroes. We first show how a met- 
ric in EF form, with u = u(r), may be put into conformal 
form. The metric Eq. J3J may be factored as follows, with 
angular dependence suppressed, 



ds 2 = 



dt 2 - dr 2 +u(dt + dr) 2 

1 — u 

(1 + u)(dt + dr)(dt - Y^~ dr ^ 
(1 + u)d(t + r)d(t - a), a i 



l-u 
l + u 



dr(30) 



The quantities t + r and t — a are termed conformal null 
coordinates since the line element is zero along lines of 
constant t + r or t — er; these thus represent radially mov- 
ing light rays. We transform to other null conformal co- 
ordinates by choosing any functions w(r + t) and v(tr — t), 
so that 

dw = w (r + t)d(r + t), dv = v'(a - t)d(a - t), (31) 
so the metric in terms of w, v is 



1 -+- u 

ds = ——dwdv — Hdwdv. 



(32) 



The metric function H may be considered to be a func- 
tion of r and t, or w and v. We may also relate the null 
coordinates to Lorentz-like coordinates, 

w = t + p, v = t — p, ds 2 = H(dr 2 — dp 2 ). (33) 

Thus r and p may be interpreted as time and radial co- 
ordinates provided that H is positive and has no singu- 
larities or zeros. 

Outside the Schwarzschild radius and the light shell 
the function a and a convenient choice for the functions 
w and v are the following 

/T + R f T \ 
dr = r + 2R\n( 1 , r > R, 
r — R \R J 



It's 



ju>a 2 r 



e- 2ar R 2aH (r - R) 



\-1aR 



2aR 



(34) 



where a, 7, a; = const. This transformation is chosen to 
make H s independent of time. In order that H s also be 



nonsingular and have no zeros we choose a = 1/2R, and 
to map the exterior region to the w s > 0, v s < quadrant 
we choose 7=1 and w = — 1. Then the transformation 
and metric functions are explicitly 

4i? 3 



R 

Note that w s and v s are dimensionless and H s is the 
square of a distance. 

For the interior of the black hole, r < R, and outside 
the light shell (see Fig. lb) we choose the opposite sign 
for the logarithm and 7 = u> = 1 in Ea. H34|) . to obtain 
the same expression as in Eq. l|35|l . which is thus valid 
throughout the Schwarzschild geometry. The expressions 
Ea. (|35|) are similar to the standard ones used to trans- 
form between Schwarzschild coordinates and KS coordi- 
nates, but differ in important ways. Lines of r =const. 
map to hyperbolae with w s v s = const.; in particular 
r = corresponds to w s v s = 1 and r = R corresponds 
to w s v s — 0. The lines t — 00 and r — R both map to 
v s — 0, while t = — 00 maps to w s = 0. (Unlike the case 
with Schwarzschild coordinates lines of constant t do not 
map to rays of v s /w s — const.!) See Fig. 0] 




FIG. 4: Collapse of a thin light shell in KS coordinates, to be 
compared with Fig. in EF coordinates. Only the spacctime 
region to the right of the line r = has physical meaning. Compare 
the light rays A B C to those in Fig. |3] 

There are many ways to transform from EF to KS 
coordinates for the Minkowski geometry inside the light 
shell; of course Minkowski geometry in EF coordinates 
is already in KS form, but not one that is useful to us. 
We must choose a transformation that joins continuously 
with the transformation Eq. I|35|l along the boundary line 
r + t = 0, and we also demand that the function H be 
continuous along that line. For that space the metric 
function u — 0, so that from Eq. 1301) a — r. Thus from 
Ea.l)32[l the metric function H is 



H m = 



-1 



w' m {r + t)v' m (r-tY 



(36) 



Equating this with H s from Eq. I|35|) along r + t = wc 
obtain a differential equation for v m 



v' m (2r) 



-1 



4R3<(0) 



re 



r/R 



(37) 
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That is, in terms of the argument, denoted x, 

v' m (x) = (const. )xe x/2R 1 

v m (x) = (const.) (a; - 2R)e x/2R . (38) 



singularity, 



For continuity we choose w r 
Schwarzschild region Eq.JJS 
stants we arrive at 



to be the same as in the 
so with appropriate con- 



Wrn = e^. 



. r — t r-t 
Vm = ( ~ ~2R^ 2R ' 

H m = 4R 3 ( — )e~^. 
r — t 



(39) 



These functions are obviously equal to those in Ea.l|35|) 
along the line r + t = 0, as desired. Note that the price 
we pay for having H s in the Schwarzschild region inde- 
pendent of time is that H m in the Minkowski region is 
dependent on time. 

The nature of the transformation Ea. (|39[l is best seen 
in terms of the mapping of some lines and points: t = oo 
maps to v = 0; t = — oo maps to w — 0; lines of constant 
r + 1 map to constant w lines; lines of constant r — t map 
to constant v lines; in particular r + 1 = maps to w = 1 
and r = t = maps to w = v = 1. Finally the origin 
r = maps to 



J/2R v 



(i 



t 

2R 



-t/2R 



-(1 + \nw m ). 



(40) 



(41) 



Henceforth we drop the subscripts on w and v. 

Figure 0] show the complete collapse process, with 
the Schwarzschild and Minkowski geometries stitched to- 
gether along the line r + 1 = 0. Only the region to the 
right of the line corresponding to r = has physical 
meaning. It is evident from the figure that the line v = 
defines a horizon, and it is also apparent that the sin- 
gularity at r — does not behave like a time indepen- 
dent spatial position. Note also that p serves as a radial 
marker, even though regions of the spacetime have neg- 
ative p. 

Figure 0] may be compared with the "pure" or eternal 
Schwarzschild geometry discussed in many textbooks 0, 
4] and shown in KS coordinates in Fig. |SJ this shows 
the "maximum analytic extension" of the Schwarzschild 
solution. The region w < 0, v < is interpreted as a 
white hole, and is absent in figure 4; also absent is the 
region w < 0, v > 0, interpreted as the "other side of the 
wormhole," and part of the interior region w > 0, v < 0. 
This illustrates that if the formation of the black hole 
by gravitational collapse is taken into account then the 
much-discussed white hole and wormhole regions are not 
present; there is no reason to expect that such regions 
occur in nature, as emphasized by Wheeler and many 
others [l7|. 




FIG. 5: The pure or eternal Schwarzschild geometry in KS coordi- 
nates. The entire spacetime region shown is given physical meaning 
in terms of a white hole region and a second exterior Schwarzschild 
region. Compare to Fig.|3] See for example references [3,4]. 



V. UNIFORM FLUID SPHERES 

The spherical shells of light used in the preceding sec- 
tions provide very simple models of gravitational collapse 
but are rather artificial and do not approximate anything 
we expect to find in nature. We now turn to a more re- 
alistic system, a sphere of uniform perfect fluid with a 
linear equation of state. This includes our app roach to 
the special case of a dust ball, the original 15] and still 
a favored system for collapse studies 0, IE W- We be- 
lieve our approach is the simplest available, because the 
Painleve-Gullstrand (PG) coordinate system is remark- 
ably well suited to the task [2(3, • 

Some of the mathematical techniques we used for light 
shells, such as use of concentric layers and the handling 
of the boundary conditions, will also prove useful for fluid 
spheres. 

The fluid collapse involves only two spacetime regions, 
the fluid interior, and the Schwarzschild exterior. A cru- 
cial step is to find a metric form (the PG form) which 
describes both regions simply, and in which the motion 
of the fluid is simple. 

The spacetime geometry corresponding to a uniform 
fluid is well known from cosmology; it is described by the 
Friedmann-Robertson- Walker metric in co-moving coor- 
dinates, which covers all of spacetime from the big bang 
onwards [3,4]. For simplicity we consider the spatially 
flat case, that is k — 0, for which the metric is 



ds 1 



dt 2 - a{t) 2 {dr 2 c + r 2 dfl 2 ). 



(42) 



(This corresponds to zero energy collapsing system.) The 
co-moving radial coordinate r c is dimensionless, while the 
scale function a(t) is a solution of the cosmological equa- 
tions, with the dimension of a length. To describe a finite 
sphere we truncate the radial coordinate at some value. 
If the fluid has a linear equation of state, with pressure 
and density obeying p — ap, then the scale function is a 
power of t, specifically 



a{t) = At 



A = const. 



(43) 
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In particular dust (or "cold matter" ) has negligible pres- 
sure so a = and n = 2/3, while radiation (or "hot 
matter") has a = 1/3 and n = 1/2. This is the range of 
normal matter. In the cosmological scenario time runs 
from the big bang at t = to infinity, but in the col- 
lapse scenario it will run from negative infinity to t = 0. 
That is, a very large fluid sphere in the far distant past 
collapses toward zero size at t = 0. 

To obtain the desired form for the metric inside the 
fluid we introduce a new radial coordinate, 

r = a{t)r c , (44) 

which is not co-moving, to obtain the metric 

ds 2 = [1 - {ra/af]dt 2 + 2{ra/a)drdt - dr 2 - r 2 dtt 2 
= [1 - (nr/t) 2 ]dt 2 + 2(nr/t)drdt - dr 2 - r 2 dtt 2 . 

(45) 

This contains the single metric function nr/t and is dis- 
tinctive in having a cross term and g rr = — 1. It has an 
infinite redshift surface where <?oo = 0, or nr/t = ±1. 
The minus sign is appropriate since we will deal with 
negative times. 

The empty region exterior to the fluid is described by 
the Schwarzschild metric in Eq.Q), which we now write 
in the form 



ds 2 = (1 - ip 2 )dt 2 s 



dr 2 



l-V' 2 



2 dn 2 , ip = ±\ 



(46) 



To make this compatible with the metric in Ea. (|45l) in 
the fluid we choose a new time coordinate that makes 
g rr = — 1. Taking 

t. = t + g(r), (47) 

we find that 

ds 2 = (1 - ip 2 )dt 2 ± 2^drdt - dr 2 - r 2 dfl 2 , (48) 
provided that g obeys 
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±- 



(49) 



As expected this transformation involves an infinite time 
stretch at the Schwarzschild radius, where ip — 1. For 
the Schwarzschild region the solution to Eg. (|49l) is 



(50) 



The metric Ea. l|48fl and the transformation Ea. (l50|l are 
those obtained by Painleve and Gullstrand [2(], l2l|. Both 
regions of spacetime are now described by the metric form 
Ea. (|48Jl . with ip allowed to be a function of both r and t; 
specifically 




fluid interior, 
Schwarzschild exterior. 



(51) 



We choose the positive signs in Eq. (|49|l and negative sign 
in Ea. H5U[) to correspond to collapse during negative time, 
and refer to the metric form and coordinates as general- 
ized Painleve-Gullstrand or simply PG. 

The PG metric has a remarkable property that is cru- 
cial to our analysis. The geodesic equations for radial 
motion of a particle in the metric ca l|48(l lead, after some 
algebra, to 



d 2 t 



/ _ dtp 



dr 



1 = (W 2 )A 2 + 2^-(^. (52) 
ds ds ds ds 



One obvious solution to the first is 



dt/ds =1, t = s — sq. 



(53) 



Thus coordinate time and proper time intervals are equal 
for such a freely falling particle, and this holds for any 
metric function ip. The second equation in Ea. (|52|l now 
becomes quite simple 



dr 
ds 



dr 
dt 



(54) 



In the Schwarzschild region this is the same as the classi- 
cal Newtonian equation for a radially falling test particle 
of zero energy. 

We will soon need the Einstein tensor and the stress- 
energy tensor. From the PG metric form it is straight- 
forward and only slightly tedious to calculate these. The 
nonzero components of the Einstein tensor are 



„ _ 2V^' i> 2 , _ 2^ 
, _ 2^' _ 4> 2 2^ 

1 _ ~ ~^2~ ~'' 



G2 
9 — \J q — 



^'_^"_^'2. (55) 



As before we define the stress-energy tensor via the field 
equations G a ^ = SnGT 1 ^; in particular the energy den- 
sity is 



J- r 



1 .2ipip' V 



8ttG r r* 
0, exterior, 
interior. 



3ra z 

STTGt' 2 



(56) 



Thus the density is uniform, as expected. 

We now join the spacetime regions inside and outside 
the fluid by demanding that the metric function ip in 
Ea. I|51fl be continuous across the fluid surface boundary. 
This gives 



r 3/2 



t = 0. 



(57) 
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The geodesic equation for a zero energy falling particle 
in the Schwarzschild exterior region is 



„3/2 



(t-t ) = 0. 



(58) 



This agrees with the motion of the surface Eq.j57|) only 
for the case of n = 2/3, that is dust. Thus a freely falling 
particle may hover at the falling surface of a dust ball, 
as it should. However, if the pressure is nonzero and 
n < 2/3 then the surface will fall more rapidly. This is 
due the stress in the surface layer, as we will see. FigureEl 
shows the collapse scenario in PG coordinates. 




FIG. 6: Collapse of a uniform fluid sphere to form a black hole in 
PG coordinates. 

It is straightforward to calculate the stress-energy ten- 
sor for the fluid surface with the same technique we used 
in section 3. In the vicinity of the surface the metric 
function tp may be written 



t K 



«) 



^e( r 3/2 + v# f) _ (59) 



This leads to a singular stress-energy tensor 

T\ = T\ = ^Lj(l - *£)6 (r - (VRt/nf/ 3 ) . (60) 

For a zero pressure fluid, a = and n = 2/3, this vanishes 
as expected. If the pressure is not zero it represents a 
surface tension, as we will now show. The surface tension 
of a fluid sphere is related to its radius and pressure by 
t = pr/2; in the present case this implies from Eq. lPt3"|) 
and Ea. ljSBT) . 

pr apr nr 3n 

T = Y- — -8^^-Y^ (61) 

which agrees with Ea. H60JI . Thus the stress-energy tensor 
indeed represents a surface tension that balances the in- 
ternal pressure of the fluid to keep it stable. This is why 
the surface falls faster than a freely falling particle. 

In summary the collapse of a uniform fluid sphere is de- 
scribed by the metric function in Ea. (|48|l and Ea. (|51|l . the 
density in Ea. (|56l) . and the surface tension in Ea. H60fl . Its 
collapse is qualitatively similar to that of the light shells, 
but rather less artificial. Our use of surface tension to sta- 
bilize the surface is a conceptually simple substitute for 



the gradual pressure gradient of a more realistic model; 
such tangential pressures are also mentioned by Singh in 
Ref. [26] and in the references contained therein. 



VI. ZERO PRESSURE FLUID SPHERES 

We next consider the collapse of a sphere filled with 
zero pressure perfect fluid, that is a dust ball. This 
problem is discussed by Landau and Lifshitz 5] using 
Lemaitre-like coordinates |34j . For simplicity and phys- 
ical clarity we use the same technique as in section 3 
for constructing thick light shells; that is we build the 
dust ball layer by layer from a sequence of thin dust 
shells. The layers are prohibited from crossing to pre- 
vent superficial singularities due to the consequent in- 
finite density [2^| . PG coordinates will again be used, 
with the metric form Ea. (|48(l containing a single metric 
function ip(r, t). The Einstein tensor for this is given in 
Ea. (|55(l . and we recall that in PG coordinates we may 
take proper time and coordinate time intervals to be 
equal along zero energy particle geodesies. 




dust shell, 
r" ! + (3-tni 2)(.0 



FIG. 7: A thin shell of dust falls into a black hole to form a larger 
black hole in PG coordinates. This is the analog of Fig.^for light. 



In analogy with section 3 we begin with a thin shell of 
dust falling onto a black hole, shown in Fig. [7| The shell 
is assumed to be very light, so that m and m' differ by a 
small Am, and the Schwarzschild radii differ by 2Am. As 
in Section 5 the metric function in the two Schwarzschild 
regions is 



- \J R/r, initial Schwarzschild region, 
-\J R' /r, final Schwarzschild region. 



(62) 



The equation of the thin dust shell is that of a geodesic 
for a zero energy particle in the Schwarzschild geometry, 
given in Eq.JSSJ). It is important that coordinate and 
proper time intervals are equal along the geodesic, so 
that the equation for the boundary is a relation between 
the coordinates r and t. 

To obtain the stress-energy tensor for the dust shell we 
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write the metric in the vicinity of the boundary as 

1>= - 




(63) 



where + = 1 as before. From the Einstein tensor in 
eq(refPGG this leads, as in section 3, to a singular energy 
density for the dust shell, 



-L n 



3 VI Vl? 7 - Vr , 

8ttG ^/2" 



5{r 3 / 2 + ^VRt). (64) 



From this we may calculate the mass of the shell by going 
to large negative times when the shell is in nearly flat 
space. 



M s = / 4nr 2 pdr = -^-(VW - VR) 



Am 



(65) 



We thus verify that the mass of the initial black hole plus 
the shell mass is equal to the mass of the final black hole. 




A/no An7 1 Am 2 ^ m f 

(a) 



(b) 



FIG. 8: In (a) a discrete sequence of dust shells forms a black hole, 
and in (b) a continuous version of the same process forms a black 
hole. These are analogs of Fig. [5]for light. 

This elementary result allows us to construct a rather 
general dust ball with non-uniform density from layers of 
thin shells. A discrete sequence of thin shells is shown in 
Fig. The initial system is a black hole of vanishingly 
small mass m with a thin shell of mass Ami collapsing 
onto it at t — to give a black hole of mass mi — mo + 
Ami, followed by more shells, and ending with a final 
shell impacting the origin at t = A to give a final black 
hole of mass m/ and Schwarzschild radius R — 2m f. In 
terms of the intermediate masses mj the metric function 
in the regions between shells is 



= —^2rnj/r, j-th Schwarzschild region. (66) 

The infinite redshift surface where ip = — 1 and goo = 
is the zigzag line in the figure. 

A continuous version of the dust shell sequence is 
shown in Figure [S^. The shells are labelled by a con- 
tinuous variable A ranging from to 1, with the total 



geometric mass inside a shell denoted by m(A); m(A) is 
the continuum analog of mj, and the metric function in 
the fluid region is thus 



ip = — y/2m(X)/r, fluid region. (67) 

The infinite redshift surface is where iji = — 1, or 

r = 2m(A), infinite redshift surface. (68) 

As in section 3 we label the shells by the total mass inside 
the shell, 



m(A) = Am/, 



(69) 



where m/ is the final mass. The time at which the A shell 
impacts the center may be chosen rather arbitrarily, and 
we denote it by h(X). For simplicity and to avoid density 
singularities it is important that the shells do not cross 
each other. This will be true if h(X) increases monoton- 
ically; that is, outer shells impact at times later than 
inner shells. Thus we choose h(X) to increase monotoni- 
cally from to A as A runs from to 1. The equation for 
the A shell is then the free fall equation for a particle of 
zero energy, with the particle reaching the origin at h(X), 
or 



h(X)] = 0. 



(70) 



As we will see, the function h(X) determines the dust ball 
density. 

If the function h(X) is specified we can use Ea. l|68[l and 
Ea. (l69|) and Eq. (|70|) to get a parametric expression for 
the infinite redshift surface; the radius and time are given 
in terms of A by 



r = RX, t = h(X) - -RX. 



(71) 



In particular for the central shell with A = 0, which im- 
pacts the center at t — 0, and the outer shell with A = 1, 
which impacts at t = A, we have for the infinite redshift 
surface 



r = t = 0, 

r = R, t = A 



central shell, 
\R, outside shell. 



(72) 



Thus the surface will move outward with time if A > 
2i?/3 - that is if the mass impacts the origin sufficiently 
slowly. Curiously the same expression occurred for the 
light shell collapse in Eq.JSHJ- 

For the special case of zero total impact time, A = 0, 
all of the shells impact the center at t = 0. This is corre- 
sponds to the uniform density fluid sphere in the previous 
section, so h(X) — clearly implies uniform density. We 
can explicitly see the relation between h(X) and the den- 
sity by solving the shell equation Eq. l|70|l for the mass 
m(A) inside the shell A at a large negative time t = —T, 
when T » h(X), to get 



m(A) 



9(T + h) 2 



9T 2 



(73) 
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Thus in the distant past the mass is proportional to r 3 , 
meaning that the dust is asymptotically uniform. As 
time progresses to smaller negative values the density 
may deviate more and more from uniform, and may be 
quite non-uniform at t = 0. 

To obtain the metric and other properties of a collaps- 
ing dust ball we first specify a function h. Then the shell 
equation Ea. H7U|) gives the shell parameter as a function 
of position, that is A = X(r,t). This is consistent because 
we allow only one shell to pass through a given point. 
With the parameter known as a function of position the 
metric in the fluid is given by Eg. 1(67(1 . With this expres- 
sion for the shell parameter the shell relation Ea. 1)70(1 is 
an identity in r, t. 

A number of interesting properties of the dust region 
may be expressed in terms of the function h(X). (Later 
we will consider a specific h(X) to illustrate further.) 
From the expression Eg. ((67(1 for the metric function ip 
and A(r, t) we may calculate the Einstein tensor from 
Eq.lJSSj, obtaining 

G °-- — > G o-^T> A-—, A-—. (74) 

For the other components we need to relate the deriva- 
tives A and A' to each other, which may be done using 
the shell relation Eq. JTUJ) . Differentiating Eq. lf7U|) we find 



dVX 
dr 

dVX 



y/r/R 



h-t + 2Xdh/dX' 

Vx 




(75) 



Thus the ratio of derivatives is independent of the func- 
tion h. With the use of Eq. 175|) it is only slightly tedious 
to show that, for any h, 



G\ = G 2 



G\ 



0. 



(76) 



Thus the pressure terms of the stress-energy tensor are 
zero, as we expect for dust. Moreover the energy density 
of the dust is 



n - T° - 

P — 1 o — 



RX' 
8irGr 2 ' 



dust region. 



(77) 



There are two specific examples of the function h(X) 
for which Eq. (|55ll is easily solved. The first example is 



h(X) = A/VA. 



(79) 



We leave it as an exercise to the reader to show that the 
metric function in the dust is then 



2A 



t+Jt 2 + (8A/3Vi?)r 3 / 2 



(80) 



In the limit of A — * this gives the same result as Ea. (|51l 
for a uniform dust ball. The energy density is 



i? 



87rGAr 3 / 2 
l/^Gt 2 ), 



t+ \/i 2 + (8A/3/R)r 3 / 2 



yjt 2 + (8A/3\/i?)r 3 / 2 
for large t. 



(81) 



This is a well-behaved positive function with no singu- 
larities or zeros for negative t. 
The second example is 



h(X) = -D/VX, 



(82) 



where D is a constant. This is singular and thus quite 
different from the first example, and is not represented 
by Fig. 8b: it is singular for A — > 0, meaning that the 
central dust layer collapsed to the center in the infinite 
past and the outer layer reached the center at time —D. 
We again leave it as an exercise to show that the metric 
function is 



if> = (2r/3i) + y/R/r(A/t), 
and that the density is 

p = D\fR/(AT:Gt 2 r z/2 ). 



(83) 



(84) 



The singularity at r = is mild in the sense that the 
mass inside a sphere goes like r 3 / 2 . As t — > the density 
approaches infinity as expected. Thus this example rep- 
resents a dust ball with a mildly singular central density 
in the distant past, which grows stronger with time until 
complete collapse, and represents a rather realistic and 
amusing system j35j. 

Finally, we briefly note some properties of light rays in 
the collapsing dust ball. From the PG metric Eq. 1)48(1 we 
may write the coordinate velocity of light as 



From this we may verify that the total mass of the dust 
ball is the mass of the final black hole, which we do by 
again going to large negative times when the dust is in 
nearly flat space, 



Mf= J^ 2 pdr=^Jx'dr = ^L. 



(78) 



Finally the stress energy tensor on the dust ball surface 
is zero according to Ea. ((61(l with n = 2/3. 



dr 
dt 



1+ip, outgoing, 
— 1 + tp, inf ailing. 



(85) 



We generally expect that tp will be zero at r = 0, except 
near t — 0, so v c = ±1 there; note that this is not true 
for the second example considered above. Also, v c — 
±1 for asymptotically large distances. On the infinite 
redshift surface ip = — 1 so that v c — for outgoing light 
and v c — — 2 for infalling light. For the center in the 
Schwarzschild region v c — — oo. This allows us to make a 
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rough qualitative sketch of some light rays in Fig. 9. (For 
A = Ea. l|85[) may be solved exactly.) Note that, as 
in the case of the thick light shell in Section 3, there is 
a last light ray out, whose path defines a global horizon 
inside of which neither particles nor photons may escape 
to infinity. 




9-2004 



FIG. 9: Some "outgoing" light rays in the collapsing dust 
ball. No ray can escape from the center after the last ray 
out (B). Compare to Fig. ThickLightShellEFHorizon for light 
shell collapse. 



VII. SUMMARY AND FURTHER STUDY 

In this paper we have tried to present a simple intro- 
duction to the dynamics of gravitational collapse, which 
we hope can provide a bridge between basic textbook 
general relativity and numerous research topics of cur- 
rent interest. The light shells in EF coordinates and the 
fluid spheres in PG coordinates provide the simplest way 
we know of to study collapse. Some simple but amusing 
extensions of these models can be made, as already men- 
tioned: one may reverse the time and study the emission 



of light and matter from white holes or other spherically 
symmetric objects, and it is easy to add a cosmological 
constant in the fluid collapse. 

The present research literature is dense with studies 
involving collapse, ranging from fundamental theory to 
astrophysical applications. A sampling of recent topics 
from the Los Alamos archives xxx.lanl.gov includes the 
following, which we paraphrase: 

• Entropy in collapse to a black hole - where does the 
information go? 

• Unitarity - is collapse consistent with unitary quan- 
tum evolution? 

• Cosmic censorship - are all singularities surrounded 
by a horizon? 

• Collapse in context of string theory, anti de Sitter 
space 

• Quantization and entropy of the surface of a black 
hole. 

• Collapse in diverse dimensions. 

• Scalar tensor (or other) gravity theories and collapse 

• Collapse with a cosmological constant included. 

• Particle production in collapse. 

• Role of pressure (radial or tangential) in collapse. 

• Stability of stars against collapse due to rotation. 

• Gravitational radiation from (non-spherical) col- 
lapse. 
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